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Abstract 

O ' 

' More than 30 miUion of high-energy muons collected with the MACRO detector at the underground 

Gran Sasso Laboratory have been used to search for flux variations of different natures. Two kinds of 
studies were carried out: search for periodic variations and for the occurrence of clusters of events. Different 
analysis methods, including Lomb-Scargle spectral analysis and Scan Test statistics have been applied to 
J> ' the data. 
I> 

oo 

^ ■ 

1 Introduction 

in 

. The high energy muon events collected by the MACRO apparatus at the average depth of 3600 m.w.e. 
represent one of the most extensive records of such kind of data. The series of these high-energy muons can 
^ , be to search for time variations of periodic and of stochastic characters, as it was done extensively by using 
arrival times of EAS [1]. These variations in the underground muon flux may be due to different causes of 
galactic, solar and terrestrial origin. The common problem for this type of searches is to determine whether 
an observed effect has occurred by chance of if it signals a departure from a pure random muon arrival. 
^ . MACRO was a multipurpose modular apparatus with 6 supermodules with scintillator detectors, limited 

streamer tubes and nuclear track detectors [2] , and studied atmospheric neutrinos [3] , aspects of CR physics 
and astrophysics [4], searched for GUT Magnetic Monopoles and other exotica [5]. Some interruptions 
' of different kinds occurred during data taking, either randomly (e.g. power outages), or regularly (e.g. 
maintenance), so appropriate statistical methods have to be applied and particular care should be used in 
choosing periods of stationary conditions. 

In the following we discuss the results of the searches for periodic variations and for time clustering of 
muon events. 
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2 Periodicity search. Spectrum analyses 

For this analysis we considered data recorded by the streamer tube system in the time interval November 
1991-May 2000 and selected the data with the following criteria: 

- run duration longer than 1 hour; 

- streamer tube efficiencies of wires and strips larger than 90% and 70%, respectively, for each module; 

- all 6 super-modules in acquisition; 

- acquisition dead time smaller than 2.5% for the whole detector. 

The total number of runs surviving these cuts was 6920 for a total number of 3.5 • 10^ muon events. 
The Fourier amplitude spectrum analysis is a powerful technique that allows a blind search for 
regular/persistent fluctuations in time series [6]. Such a technique, however, requires the input data to 
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Figure 1: Lomb power as a function of the Logio of the frequency [days~^] for experimental data (upper 
panel). Note the high peak at ~ —2.56 corresponding 365 days (seasonal flux variation). In the lower panel 
the results of a Monte Carlo simulation having the same noise level of the real data with seasonal, solar 
diurnal and sidereal waves [9] added. 



be sampled at evenly spaced intervals; data gaps of variable length and occurring randomly in the serie 
produce spurious contributions to the power that can mimic the presence of a periodicity. The Lomb- 
Scargle method [7] has been developed to mitigate this effect even in the case of very long data scries. 
Moreover, as indicated in [8], it allows to evaluate the significance of the "peaks" (signal) with respect 
to a null hypothesis. The muon events were binned in 15 min time intervals and bins deviating by more 
than 3a from the monthly average rate were discarded. The total number of time bins used was 160242 
corresponding to 58% of the whole sample. 

The results of our analysis are shown in Fig. 1. We compare the spectrum obtained for the real data 
with a Monte Carlo simulation having the same noise level and time intervals distributed according to the 
sequence of the original series. The seasonal, solar diurnal and sidereal waves [9] were also added in the 
serie. The spectrum of real data shows a power distribution similar to what observed in other cosmic ray 
data series [6], i.e. a low frequency spectrum whose power decreases with frequency"^. The most striking 
feature of the spectrum is the large peak at ~ —2.56 corresponding to the seasonal flux variation. Fig. 2 
shows a frequency region around the solar diurnal frequency where we have also indicated the frequencies 
corresponding to the sidereal and anti-sidereal waves. The peak at frequency 1 day~^ has a statistical 
significance of ~ 2.3(t; the statistical significance assuming an oscillatory behaviour is ~ 3Aa. Note that 
peaks corresponding to the diurnal and sidereal variations are observed, but peaks of similar size (or even 
larger) are present elsewhere in the spectrum. The claim that the sidereal and solar diurnal waves are real 
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Figure 2: The frequency region around the solar diurnal wave. The arrows mark its position and the sidereal 
and anti-sidereal peaks. 



is based upon its occurrence at a frequency of "a priori" interest and on the stability of its amplitude and 
phase with time. We find that the amplitudes and the probabilities for the null hypothesis computed by 
the two methods are in fair agreement with the ones obtained using a standard "folding" method [9] . 



3 Burst search: time interval distribution 

The first method used in searching for correlations in the arrival times was the study of the time interval 
distribution. For each muon arriving at time to we calculated the distribution of the time interval elapsed 
between the first muon to and the next five muons: tj-to, i=l,.--,5, see Fig. 3. A complete analysis was 
published in [10]. Here we report the results for the direction bands with 0° < RA< 360° and 25< decl 
< 50° that include the Cyg X-3 region. The experimental distributions show some deviations from the 
muon random arrival expectation. The probability computed using the Kolmogorov-Smirnov test shows 
some disagreement (prob=0.38) but the available statistics is too poor to reach clear conclusions. 

4 Burst search: scan statistics 

Scan statistics is a powerful method to search for bursts of events. It is a bin-free method and it provides 
unbiased results when data are analysed (see [11] and references therein). We used scan statistics in the 
following way: for each run i, let [Ai, Bi] be the time interval ranging from the start and the end of the run. 
We open a "time window" of fixed length w and scan the interval [j4i,_Bi] counting the number of events 
falling inside w. ki is the maximum number of events recorded during the scan. Finally, for each run. we 
compute the probability Pi that a statistical fluctuation would produce a burst of events as large as ki. The 
only choice is the size of w. We tried different sizes (w = 30 s, 5 min and 15 min) and for each of them 
we analysed the probability distribution Pi, i=l, Nj.„„. In Fig. 4 we show the probability distribution for 
the 6113 runs surviving our cuts: = 30 s above, 5 min at the centre and 15 min below. No significant 
deviations from the null hypothesis is found; we also inspected unusual runs with probabilities smaller than 

5 • 10"^ and we found that the "bumps" of events were concentrated near the beginning or end of the runs. 
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Figure 3: Time interval distributions (from top left to bottom left clockwise) for t2-to, ts-to, ts-to and t4-to, 

respectively, for a cone of arrival directions with declination 25° < S < 50°. The dashed lines represent the 
fits to the Erlangen gamma function of order 2, 3, 5, 4, respectively (see [10]). 
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Figure 4: Scan statistics probability distributions for the selected runs. In the upper panel a time window 
w = 30 s was used; w = 5 min in the central panel and w = 15 min in the bottom panel. 
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5 Conclusions 



We analysed the time series of MACRO muons using two complementary approaches: search for periodicities 
and search for bursts of events. The Lomb-Scargle method was used in the first case, scan statistics in the 
second. The two techniques complete early analyses performed with "folding" methods in searching for 
periodicities and time differences for burst events. No deviations from the expected distributions were 
found. 
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